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Abstract 

We investigate a system related to a particular isothermal gas-solid chromatography process, 
called "Pressure Swing Adsorption", with two species and instantaneous exchange kinetics. This 
system presents the particularity to have a linearly degenerate eigenvalue: this allows the velocity 
of the gaseous mixture to propagate high frequency waves. In case of smooth concentrations with 
a general isotherm, we prove L 1 stability for concentrations with respect to weak limits of the 
inlet boundary velocity. Using the Front Tracking Algorithm (FTA), we prove a similar result for 
concentrations with bounded variation (BV) under some convex assumptions on the isotherms. In 
both cases we show that high frequency oscillations with large amplitude of the inlet velocity can 
propagate without affecting the concentrations. 

Key words: systems of conservation laws, boundary conditions, BV estimates, entropy solu- 
tions, linearly degenerate helds, convex isotherms, Front Tracking Algorithm, waves interaction, 
geometric optics. 
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1 Introduction 

"Pressure Swing Adsorption (PSA) is a technology used to separate some species from a gas under 
pressure according to the molecular characteristics and affinity of the species for an adsorbent 
material. Special adsorptive materials (e.g. zeolites) are used as a molecular sieve, preferentially 
adsorbing the undesired gases at high pressure. The process then swings to low pressure to desorb 
the adsorbent material" (source: Wikipedia). 

A typical PSA system involves a cyclic process where a number of connected vessels containing 
adsorbent material undergo successive pressurization and depressurization steps in order to produce 
a continuous stream of purified product gas. We focus here on a step of the cyclic process, restricted 
to isothermal behavior. 

As in general fixed bed chromatography, each of the d species (d > 2) simultaneously exists under 
two phases, a gaseous and movable one with velocity u(t, x) and concentration a(t,x) or a solid 
(adsorbed) other with concentration qi(t,x), 1 < i < d. We assume that mass exchanges between 
the mobile and the stationary phases are infinitely fast, thus the two phases are constantly at 
composition equilibrium: the concentrations in the solid phase are given by some relations = 
9*(ci, Cd) where the functions q* are the so-called equilibrium isotherms. A theoretical study of 
a model with finite exchange kinetics was presented in [5] and a numerical approach was developed 
in 0. 
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In gas chromatography, velocity variations accompany changes in gas composition, especially in 
the case of high concentration solute: it is known as the sorption effect. In the present model, the 
sorption effect is taken into account through a constraint on the pressure (or on the density in this 
isothermal case). See [31] and [35] for a precise description of the process and [10] for a survey on 
various related models. 

The system for two species (d = 2) with three unknowns (u, ci, C2) is: 

d t (ci+q* 1 (c ll c 2 )) + d x (uc 1 ) = 0, (1) 
d t {c2 + q* 2 (c 1 ,C2)) + d x (uc 2 ) = 0, (2) 
ci+c 2 = p(t), (3) 

with suitable initial and boundary data. The function p represents the given total density of the 
mixture. The experimental device is realized so that it is a given function depending only upon 
time and in the sequel we assume that p = 1 (which is not really restrictive from a theoretical 
point of view). First existence results of large solutions satisfying some entropy criterium in the 
case of two chemical species were obtained in [5] [5] • In the previous system, it appears that we can 
expect strong singularities with respect to time for the velocity u. For instance, let c\(t, x) = c 1 be 
a constant, c 2 (t,x) = 1 — c 1; u(t,x) = u b (t) where u b is any L°° function, then (ci,C2,u) is a weak 
solution of (|TJ) , d2J) , ((3]) . So we can build solutions with a strong oscillating velocity for this system. 
Furthermore high oscillations of the incoming velocity u b slightly perturb the concentration as we 
will see. Notice that we seek positive solutions (01,02), thus, in view of ([3]) with p = 1, ci, c% must 
satisfy 

< ci, c 2 < 1. 

We use the following notations, introduced in [3]: we set c = c\ € [0, 1] and 



Qi( c ) = <7*( c , !- c ), i = l,2, 
h(c) = q 1 (c)+q 2 (c), 
1(c) = c + qi(c). 

Adding (TIJ and © we get, thanks to ([3]): 

dt(qi(c) + (72(c)) + d x u = 0, 

thus our purpose is to study the following system: 

r d t i(c) + d x (uc) = 0, 

[ d t h(c) + d x u = 0, 

supplemented by initial and boundary data: 

c(0,x) = c (x) G [0, 1], x > 0, 
< c(t,0) = c b (t)e[0, 1], t>0, (5) 
u(t,0) = u b (t) > 0, t > 0. 

Notice that we assume in §5§ an incoming flux at the boundary, i.e. Vt > 0, Ufc(i) > 0. In the case 
where the first species is inert, that is q\ = 0, the I function reduces to identity. 

System has a null eigenvalue as the system exposed in [4 \. but, instead [4], we cannot 
reduce this system to a single equation for general solutions with shocks. In [3] is studied another 
interesting 2x2 system with a linearly degenerate eigenvalue which modelises some traffic flow. As 
in [17l HH [TB] EH [30] > the zero eigenvalue makes possible the existence of stratified solutions or the 
propagation of large-amplitude high frequency waves. Usually, for genuinely nonlinear conservation 
laws, only high oscillating solutions with small amplitude can propagate: see for instance [22l I14j . 

In this paper we prove, for large data, that the velocity is a stratified solution in the following 
sense: u(t, x) = u b (t) v(t, x) where v is as regular as the concentration c and more than the boundary 
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data tit,. This decomposition for the velocity allows high oscillations with large amplitude for 
velocity to propagate, without affecting the concentration. For this quasilinear system we have 
propagation of high oscillations with large amplitude for velocity as in a semilinear system, see for 
instance [24, 25j[27], and we have strong profile for u with double scale as in [26] , 

This also permits to pass to the weak limit for u at the boundary and to the strong limit in the 
interior for the concentration. For the smooth case, we have no restriction on the isotherms, but 
for the realistic case with shock-waves, we restrict ourselves to the classical treatment of hyperbolic 
systems: eigenvalues are linearly degenerate or genuinely nonlinear. Furthermore we obtain better 
interaction estimates when the shock and rarefaction curves are monotonous. It is the case for 
instance for an inert gas and an active gas with the Langmuir isotherm. We conjecture that our 
result is still valid for general isotherms with piecewise genuinely nonlinear eigenvalue. 

The paper is organized as follows. In Section [5] we recall some basics results from [5] concerning 
hyperbolicity, entropies, weak entropy solutions of System Q. 

In Section [3l we study the case where concentration is smooth and the velocity is only 

In the remainder of the paper we study the case with only BV concentrations. In short section 
[4] we briefly expose the Front Tracking Algorithm (FTA) for System Q. 

Sectionals devoted to the study of both shock and rarefaction curves. We state the assumptions 
that we need to perform estimates with the Front Tracking Algorithm. These assumptions restrict 
us to convex (or concave) isotherms and we give some examples from chemistry. We obtain the 
fundamental interaction estimates in Section [6] and BV estimates for v in Section [7] Finally, we 
obtain strong stability for concentration with respect to weak limit on the boundary velocity in 
Section [S] 



2 Hyperbolicity and entropies 



In order the paper to be self contain, we recall without any proof some results exposed in [9] . 

It is well known that it is possible to analyze the system of Chromatography, and thus System 
([3]), in terms of hyperbolic system of P.D.E. provided we exchange the time and space variables and 

u 



u > 0: see [32] and also [34] for instance. In this framework the vector state will be U 

m 

where m = u c is the flow rate of the first species. In this vector state, u must be understood as 
u p, that is the total flow rate. 

In the sequel, we will make use of the function / = q\ ci — qi Ci introduced by Douglas and al. 
in [28], written here under the form 

/(c) = qi c 2 - q 2 ci = q\{c) - ch(c). (6) 

Any equilibrium isotherm related to a given species is always increasing with respect to the corre- 

dq* 

sponding concentration (see 281) i.e. — - > 0. Since c — c\ and C2 = 1 — c, it follows: 

aci 

q[>0> q' 2 - (?) 

Let us define the function H by 

H{c) = l + (l-c)q' 1 (c)-cq' 2 (c) = 1 + q[{c) - ch'{c). (8) 
From ([7]), H satisfies H > 1 and we have the following relation between /, H and h: 

/"(c) = H'(c)-h'(c). 
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2.1 Hyperbolicity 

Concerning hyperbolicity, we refer to [20l [36], [37] . System ((4]) takes the form 

h(m/u) 

I(m/u) 
The eigenvalues are: 



d x U + <9 t $(t/) = with U = [ U m ) and $([/) = ( ) . (9) 



and A 



H(c) 



thus in view of |8j the system is strictly hyperbolic. The zero eigenvalue is of course linearly 

degenerate, moreover the right eigenvector r = ( , , ?l x ] associated to A satisfies dX ■ r = 

V 1 + 9i( c ) / 



,2 



/ (c), thus A is genuinely nonlinear in each domain where /" =/= 0. 



Proposition 2.1 (|9j Riemann invariants) 

System admits the two Riemann invariants: 

c and w = ln(u) + g(c) = L + ,g(c), where g'(c) — JjfTY an ^ ^ = l n ( w )- 
Furthermore this system can be rewritten for smooth solutions as: 

d x c + ^ d t c = 0, d x (\n(u) + g(c)) = d x w = 0. (10) 

2.2 Entropies 

Dealing with entropies, it is more convenient, as shown in [9], to work with the functions 

G(c) = exp(p(c)), = exp(w) = u G(c). 
Notice that G is a positive solution of HG' + h'G = 0. 

Denote E(c, u) any smooth entropy and Q = Q(c, u) any associated entropy flux. Then, for smooth 
solutions, d x E + dtQ — 0. Moreover: 

Proposition 2.2 (|9j Representation of all smooth entropies) 

The smooth entropy functions for System are given by 

E(c,u) — 4>{w) + utj}(c) 

where (j) and ip are any smooth real functions. The corresponding entropy fluxes satisfy 

Q'(c) = h'{c)^{c)+H{c)i)'{c). 

Moreover, in [7j, the authors looked for convex entropies for System (i.e. System (fj| written 
in the (it, m) variables) in order to get a kinetic formulation. The next proposition gives us a family 
of degenerate convex entropies independently of convexity of / or of the isotherms. 

Proposition 2.3 (|9j Existence of degenerate convex entropies) 

If if) is convex or degenerate convex, i.e. tp" > 0, then E = uip(c) is a degenerate convex entropy. 

There are some few cases (water vapor or ammonia for instance) where the isotherm is convex. 
There is also the important case with an inert carrier gas and an active gas with a concave or 
convex isotherm (see [5J |H1 EE])- I n these cases, the next proposition ensures the existence of A- 
Riemann invariants which are also strictly convex entropies. In such cases, w is monotonous with 
respect to x for any entropy solution. 
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Proposition 2.4 (|9j When A-Riemann invariant is a convex entropy) 

There are strictly convex entropy of the form E = <j)(w) if and only if G" does not vanish. 

More precisely, for a > 0, E a (c, u) — u a G a (c) is an increasing entropy with respect to the Riemann 

invariant W . It is strictly convex for a > 1 if G" > and for a < 1 if G" < 0. 

Unfortunately, when G has an inflexion point such system does not admit any strictly convex 
entropy. When one gas is inert, it is always the case if the sign of the second derivative of the 
isotherm changes. See for instance [S] for the BET isotherm. 

Remark 2.1 In general, System Op is not in the Temple class. It is the case if and only if f" 
does not vanish and d x W = for all entropy solution (fUj/). For instance, System (CJ) with two 
linear isotherms is in the Temple class. 

Proposition 2.5 (|9j Non Existence of strictly convex entropy) 

If sign of G" changes then System ffl) does not admit strictly convex smooth entropy. 



2.3 Definition of weak entropy solution 

We have seen that there are two families of entropies: utp(c) and <j>(uG(c)). 

The first family is degenerate convex (in variables (u, uc)) provided ip" > 0. So, we seek after weak 

entropy solutions which satisfy d x (uip(c)) + d t Q(c) < in the distribution sense. 

The second family is not always convex. There are only two interesting cases, namely ±G"(c) > 

for all c <E [0, 1]. When G" > and a > 1, we expect to have d x (u G(c)) a < from Proposition ^. 41 

But, the mapping W i— > W a is increasing on R + . So, the last inequality reduces to d x (uG(c)) < 0. 

In the same way, if G" < 0, we get d x (u G(c)) > 0. 

Now, we can state a mathematical definition of weak entropy solutions. 

Definition 2.1 Let be T > 0, X > 0, u e L°°((0,T) x (0,X),R+), < c(t,x) < p = 1 for almost 
all (t,x) G (0, T) x (0, X). Then (c,u) is a weak entropy solution of System OJ)-f3P with respect 
to the family of entropies uip(c) if, for all convex (or degenerate convex) ip: 

^( U ^( C )) + |Q( C ) < 0, (11) 
m V'([0,T[x[0,X[), where Q' = Hip' + h'tp, that is, for all 6 D([Q,T[x[0, X[): 

X t-T pT pX 

(uip(c)d x c/) + Q(c)dt(j>) dtdx + / u b (t)ip(c b (t)) cf>(t,0) dt + / Q(c (x)) (j){Q,x) dx > 0. 



Remark 2.2 If ±G" > then uip = ±uG(c) is a degenerate convex entropy, with entropy flux 
Q = 0, contained in the family of entropies uip(c). So, if G" keeps a constant sign on [0, 1], (c, u) 
has to satisfy: 

±-0-(u G(c)) < 0, if ± G" > on [0, 1] . (12) 
ox 

Notice that the entropies uip(c) and the entropy uG(c) are linear with respect to the velocity u. 



2.4 About the Riemann Problem 

The implementation of the Front Tracking Algorithm used extensively from Section |4] requires some 
results about the solvability of the following Riemann problem: 

d x u + d t h(c) — 0, 

(13) 

d x {uc)+d t I{c) = 0, 
c M =0-610,1], *>0, I %l [0 Q f' t>0. (14) 

We are classically looking for a selfsimilar solution, i.e.: c(t,x) = C(z), u(t,x) — U(z) with 
z = — > 0. The answer is given by the three following results ([S]). 
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Proposition 2.6 Assume for instance that < a < c < c + < b < 1 and f" > in ]a, b[. Then 
the only smooth self-similar solution of \ 13\) - [l4\ ) is such that : 

< z < z-, 
z- < z < z + , (15) 
z+ < z, 

where z + = — — - , z~ = z + e~*( c+ ' with $(<:) = j ^ ^} <i£. Moreover u~ = — — - and U is 



C(z) = 


c_, 


dC 


H{C) 


dz 


Z f"{cy 


C(z) = 





given by: 




U-, < z < z- 

H(C(z)) 



< 

z_ < z < z + , (16) 



z 

u + z + < z. 

Proposition 2.7 If (c _ ,c + ) satisfies the following admissibility condition equivalent to the Liu 
entropy- condition ('291): 

for all c between c~ and c + , — — ^ — — - < — — -, 

c+ — c c — c 

then the Riemann problem H3\) - (L^ is solved by a shock wave defined as: 

where u~ and the speed s of the shock are obtained through 

[f] , 1 + [f] l + h+ 



where [c] = c+ - c~ , [/] = /+ - /" = /(c+) - /(c~), ^+ = /i(c+), ft- = ft(c"). 

Proposition 2.8 Two states U~ and U + are connected by a contact discontinuity if and only if 
cT = c + (with of course u~ ^= u + ), or c~ ^ c + and f is affine between c~ and c + . 

It appears from these results that we can build a weak entropy solution of the Riemann problem 
(|13[) - (fTl| in a very simple way (see [S]), similar to the scalar case with flux /, for any data. In 
particular, if /" has a constant sign (which is the framework in Section 3]), the Riemann problem 
is always solved by a simple wave. 



3 Case with smooth concentration 

System Q has the strong property that there exist weak entropy solutions with smooth concentra- 
tion c on (0,T) x (0, X) but not necessarily smooth velocity u, for some positive constants T and 
X. Furthermore, c is the solution of a scalar conservation law. 

3.1 Existence of weak entropy solutions with smooth concentration 

For this section, we refer to [16] , [15] . We have a similar result in [8] but only with smooth velocity. 
Here, we obtain by the classical method of characteristics existence and uniqueness of a weak 
entropy solution with smooth concentration and only L°° velocity. 
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Theorem 3.1 (Unique weak entropy solution with smooth concentration) 

Let T a >0,X>0,c € W^^X], [0,1]), c b G ^°°([0, T ], [0,1]), lnu b G L°°([0, T ], R). 
7/co(0) = Cft(0) i/ien i/iere exists T g]0, To] suc/i t/iai System {^p-{5|) admits a unique wea/c entropy 
solution (c, it) on [0,T] x [0,X] with 

cG ^'"([O.T] x [0,X],[0,1]), lnu £ L DO ([0 ! T],^ DO ([0,X],R)). 

Furthermore, for any ip £ C 1 ([0, 1],R), setting 

F'(c) = (H(c) G(c))' 1 and Q' =Hip' + h' ip, 

(c, u) satisfies: 

d x {uip(c))+d t Q(c)=0, d x (uG(c))=0, (18) 

a t c+« t (()G(c t (())9/(c)-0. (19) 

Proof: we build a solution using the Riemann invariants and we check that such a solution is 
an entropy solution. Next, we prove uniqueness. 

Using the Riemann invariant W = u G(c) (d x W = 0) and the boundary data we define u by: 

_ u b {t)G(c b {t)) 
Kh ' G(c(t,x)) ' 

so u is smooth with respect to x. Then, the first equation of (|10| can be rewritten as follows: 



d t c + fid x c = 0, with ^ = A = = — = c). (20) 

We solve (|20p supplemented by initial-boundary value data (cq, c&) by the standard characteristics 
method. Let us define, for a given (t,x), X(-,t,x) as the solution of: 

dX(s t, x) _ x ^ ^ ^ ^ _ ^ 

as 

dc 

Since — (s, X(s, r, x)) = from (|2H)) . we have 
as 

X(s,r, a;) = x — b(s,r) F'(c{t, x)) with 6(s,r) = / u b (a) G(c b (a)) da. 



Now, for some T G [0, To] defined later on, we split £1 = [0, T] x [0, X] according to the characteristic 
line T issuing from the corner (0, 0), i.e. we define the sets fi ± = {(<, x) G fi, ±(x — 0, 0)) > 0}. 
Since d x X(t,0,x) = 1 - b(t, 0) F"{cq[x)) d x c {x), 6(0,0) = and 6(.,0) G W 1 ' 00 ^), the mapping 
a; 1 > X(t,0,x) is a Lipschitz diffeomorphism for < t < T with T G]0,T)] small enough. Then 
we define on Q + , for each i G [0,T], £(t, x) such as X(t,0,£(t,x)) — x. Then we have c(t,x) — 
Co(£(t,x)) on Q + . Furthermore dt£, — ~d s X/d x X and thus c is Lipschitz continuous in time and 
space on 0+ . 

We work in a similar way on Cl~ and c G W Since c is continuous on T from the compa- 

tibility conditions c (0) = c b (0) we have c G W /1: °°(fi). 
By construction (c, it) satisfies (jTOJ) rewritten as follows: 

9a, lnu = — d x g{c), ud x c + Hdtc = 0. 

These equations imply: 

d x u— —ud x g(c) = —ug'(c)d x c— —ug'{c) ( — dtc ) = —h'(c)d t c— —d t h(c). 
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Now we check that (c, u) satisfies (IT51) . Let ip be a C 1 function. Using the identity Q' = /i'^> + Hip' 
and the previous equations we have: 

d x (uip(c)) + d t Q(c) — ipd x u + uip' d x c + Q' d t c = ip (d x u + h! d t c) + ip' (ud x c + H d t c) 
= V x + ip' x = 0. 

By the way (TT8|) implies (fTTj) . so (c, u) is an entropy solution of System ||3J|. 
We now prove the uniqueness of such a weak entropy solution. 

Precisely, if c G W^°°([0,T] x [0,X],[Q,1]) and lnw G L°°((0, T), W 1,o °(0, X)) satisfy HI]) in 
T)'([0, T[x [0, X[) with initial-boundary data co,c b ,u b then we show that (c, u) is necessarily the 
previous solution built by the method of characteristics. 

Choosing the convex functions ip(c) — ±1 and ip(c) — ±c we obtain ([4]). The main ingredient to 
conclude the proof is the fact that u admits a classical partial derivative only with respect to x. 
Thus classical computations with smooth functions to obtain (fT0|) as in the proof of Proposition ^. II 
are still valid. Now (c, u) satisfies (fT0|). which implies from the beginning of the proof of Theorem 
13. II that (c, u) is our previous solution. □ 



Remark 3.1 

1. Notice thatT,X are only depending on \\ m(u(,)||i,<jo } ||ci,||wi,<n, ||co||v^ 1 . < »- Thus, if q< e <i 
is a sequence of boundary velocity data such that (lnuf) is uniformly bounded in L°°(0,To), 
and if (c§), (c|) are some initial and boundary concentration data uniformly bounded in W 1,a ° 
with the compatibility condition at the corner Co(0) = Q,(0), then there exist T > and X > 
and Lipschitz bounds for c e ,lnu e on [0, T] x [0, X] independent of e. 

2. As in JSC, we have a global solution with smooth concentration if X is genuinely nonlinear (for 
instance an inert case and a Langmuir isotherm) , with monotonicity assumptions on cq and 
c b . 



3.2 Strong stability with respect to velocity 



In case of a Lipschitz continuous concentration, we now give a strong stability result for the con- 
centration with respect to a weak limit of the boundary velocity. 

Theorem 3.2 (Strong stability for smooth concentration) 

Let be T > 0, X > 0, c G W 1 '°°([0,X], [0,1]), c b G W 1 ' 00 ^, T ], [0,1]) such that c (0) = c b (0), 
and (lnu|)o<e<i a, bounded sequence in L°°(0,Tq). Then, there exists T g]0,To[ such that System 
admits a unique weak entropy solution (c e ,u e ) with c e G W /1 ' oo ([0,T] x [0, X], [0, 1]), In w e G 
L°°([0, T], W 1,o °([0, X], K)) and with initial and boundary values: 



c £ {0,x) = Co (x) G [0, 1] , x > 0, 
c E (t,0) = Cb(i) G [0, 1], £>0, 
u £ (t,0) = u e b (t)>0, t>0. 

If (u|) converges towards Tif, in L°°(0,Tq) weak—* when e goes to 0, then (c e ) converges in 
L°°([0,T] x [0, X]) towards the unique smooth solution of 

d t c + u b (t) G(c b (t)) d x F(c) = 0, c(t, 0) = c b (t), c(0, x) = c (x). 

Furthermore we have: 

G(c b (t)) 



(21) 



(22) 



lim 



u e (t, x) 



m 



G(c(t,x)) 



= 0. 



L°°([0,T]x[0,X]) 



Proof: thanks to Theorem l3.U there exists T > such that System with initial and boundary 
values (|21|) admits the unique weak entropy solution (c e ,u 6 ) with smooth concentration in the 



S 



previous sense on [0,T] x [0, X}. 

Since (c £ ) is bounded in W 1,oc , up to a subsequence, (c £ ) converges strongly in L°° to c. Using (fl~9|) 
in conservative form, we can pass to the limit and get (|22[) . Problem ([22|) has a unique solution by 
the method of characteristics. Thus, the whole sequence (c e ) converges. We recover the last limit 
for u £ thanks to d x (uG(c)) = 0. □ 

Notice that if Ub is a constant function for instance u £ b (t) — Ub(t/e) with Ub periodic we can 
compute the concentration with only using a constant velocity (the mean velocity) as in liquid 
chromatography. 

An example from geometric optics: if = w& ^t, -\ where Ub{t,6) 6 L°°((0, T), C°(K/Z)), 



inf Ub > 0, we have a similar result with Equation 
profile U: 



for c where Ub{t) 



Ub(t, 9) d9 and a 



lim 



u £ (t, x)-U 



t 

t 1 x } — 

e 



= where U(t, x, 9) = Ub(t, I 



G(c b (t)) 
G(c(t,x))- 



4 Front Tracking Algorithm 

In Section [31 where c is smooth and lnu;, is in L°°(0,T), we have seen that there exists a smooth 
function v such that 

u(t,x) = u b (t)v(t,x). (23) 

Furthermore c satisfies the scalar conservation law (fT5)) . For only BV data we cannot expect to 
obtain such a scalar conservation law for the concentration, except in the case of linear isotherms. In 
that case, the scalar conservation law (fH?|) and System (0| have the same solution for the Riemann 
Problem, but linear isotherms are of a poor interest from Chemical Engineering point of view. The 
first interesting case is the case with an inert gas and a Langmuir isotherm, first mathematically 
studied in [5]. 

Nevertheless we guess that (|2"3"]l is still true with v E BV. From we have yet obtained BV 
regularity with respect to x with a Godunov scheme. To get BV regularity with respect to t we 
will use a more precise algorithm to study wave interactions, namely a Front Tracking Algorithm 
(FTA). 

The Front Tracking method for scalar conservation laws was introduced by Dafermos, [TH]. The 
method was extended to genuinely nonlinear systems of two conservation laws by DiPerna [21j . For 
our purpose, we do not use the generalisation to genuinely nonlinear systems of any size by Bressan 
[H] or Risebro [33]. 

The FTA is much more complicated when an eigenvalue is piecewise genuinely nonlinear, see 
[21[TJ[53]. Then, we restrict ourselves to the case where A is genuinely nonlinear, which allows us 
to treat some relevant cases from the point of view of chemical engineering like an inert gas with a 
Langmuir isotherm, two active gas with a binary Langmuir isotherm for instance. For this purpose 
we work in the framework exposed in the recent and yet classical Bressan's Book [13] . In this 
framework we assume /" > 0, then a Riemann problem presents only two waves: 

1. a contact discontinuity with speed 0, 

2. a rarefaction wave with speed A > 

or a shock wave with speed between A - and A + , characteristic speeds associated to the left 
and right states, respectively. 

Let be 8 > 0. A ^-approximate Front Tracking solution of System ^ is a pair of piecewise constant 
function c s (t,x),u s (t,x), whose jumps are located along finitely many straight lines t — t a (x) in 
the t — x plane and approximately satisfy the entropy conditions. For each x > and ip" > 0, one 
should thus have an estimate of the form: 

J2([u 5 ^c s )}-^[Q(c s )])(t a ,x) < 0{8), (24) 
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where [u] = u + —u is the jump across a jump line, and the sum is taken over all jump for x fixed. 
Inequality (|24p implies that (c s ,u s ) is "almost an entropy solution": 

d x u s ^(c s ) + d t ^(c s ) <0{5). (25) 

That's enough to get an entropy solution "issued from FTA" when S goes to zero. 

Since we want to only use piecewise constant functions, it is convenient to approximate a 
continuous rarefaction wave by a piecewise constant function. For this purpose, the rarefaction 
curve is dicretized with a step of order 8 and then (|24|) still holds. 

We now briefly describe an algorithm which generates these Front Tracking approximations. 
The construction starts on the initial line x = and the boundary t — by taking a piecewise 
constant approximation of initial value Cb(t),Ub(t) and boundary values cq(x). Let t\ < ••• < t^, 
x\ < ■ ■ ■ < xm be the points where initial-boundary values are discontinuous. For each a = 
1, ••• ,iV, the Riemann problem generated by the jump of initial constant values at (t a ,x = 0) 
is approximately solved on a forward neighborhood of (t a , 0) in the t — x plane by a function 
invariant on line t — t a — ax, for all positive a, and piecewise constant. Notice that the boundary 
is characteristic, then we have only one wave associated with the speed A in the corner (0,0). 

The approximate solution (c s ,u s ) can then be prolonged until X\ > is reached, when the 
first set of interactions between two wave-fronts takes place. If x\ > x\ we first have to solve the 
characteristic boundary Riemann problem at (t = 0, x = x\). Since (c 5 ,u 5 ){.,x x ) is still a piecewise 
constant function, the corresponding Riemann problems can again be approximately solved within 
the class of piecewise constant functions. The solution is then continued up to a value X2 where 
the next characteristic boundary Riemann problem occurs or the second set of wave interactions 
takes place, and so on. 

According to this algorithm, contact discontinuity fronts travel with speed zero, shock fronts 
travel exactly with Rankine-Hugoniot speed, while rarefaction fronts travel with an approximate 
characteristic speed. However, one exception to this rule must be allowed if three or more fronts 
meet at the same point. To avoid this situation, we must change the positive speed A of one of the 
incoming shock fronts or rarefaction fronts. Of course this change of speed can be chosen arbitrarily 
small and we have again Inequality (|24|) . 

Notice that, for 2x2 system the number of wave- fronts cannot approach infinity in finite x > 0. 
DiPcrna shows in [21] that the process of regenerating the solution by solving local Riemann 
problems yields an approximating solution within the class of piecewise constant functions that is 
globally defined and that contains only a finite number of discontinuities in any compact subset of 
the t — x quarter plane t > 0, x > 0. We then do not consider non-physical fronts as in |13j for 
general n x n systems with n > 3. 



5 About the shock and rarefaction curves 

In this section we state our assumptions to use the FTA with large data. Precisely we work in 
classical hyperbolic case, namely, eigenvalues are linearly degenerate or genuinely nonlinear. We 
assume: 

A = — — is genuinely nonlinear i.e / is convex on [0, 1]. (26) 

Actually A is genuinely nonlinear for /" ^ 0, but since / = c\q2 — c^qi (see (JS])) we can assume 
that /" > exchanging the gas labels 1 and 2 if necessary. 

Our analysis of wave interactions in Section[B]is more precise with monotonous A-wave curves, then 
we also assume: 

A-wave curves are monotonous. (27) 

To state precisely this last assumption let us introduce some notations. Let (c_,L_) be a left 
constant state connected to (c+,L+) a right constant state by a A-wave curve. In the genuinely 
linear case, with Assumption (|26p , A-wave curve is a rarefaction curve with c_ < c + or a shock 
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curve with c_ > c+. The sign of [c] = c + — c_ comes from the general study of the Riemann 
problem in [9]. From the Riemann invariant w = lnu + g(c) and the Rankine-Hugoniot conditions 
a A- wave curve can be written as follows (see [9]): 

[L] = L + -L_= In u+ -ln„_ = T(c +) c_) = ( ^ = " ^ < C + (28) 



5(c+, c_) else 
We give an explicit formula for S* in Lemma |5. II 

Notice that we use only one Riemann invariant, namely c, to write A- wave curves. Indeed L = Inu 
and c have quite different behavior as seen in [El [9] and this paper. Furthermore we can give some 
simple criterion to have monotonous A- wave curves. For instance, as g' = —h'/H, the rarefaction 
curve is monotonous if and only if h is monotonous. A chemical example, investigated in [Sj, is the 

K2C2 

case of an inert gas (qi — 0) and an active gas with a Langmuir isotherm: q 2 (c 2 ) = Q2 77 — • 

1 + K 2 c 2 

For this case we have 

/">0, h'<0, |±> >|^-. (29) 

oc- ac + 

The first condition of (|2"9"]1 gives us (|2l)]) and the last one gives us (f2"T)) . 
Notice that if we exchange labels 1 and 2 for gas, Inequalities (|29p simply become: 

r <„<„. ^<o<^. 

oc- oc+ 

Let us give some isotherm examples such that ([2U]l is satisfied. 

Proposition 5.1 For the following examples, Assumptions \26\) , p?7| ) are valid: 

1. one gas is inert: q\ — 0, and the other has a concave isotherm: q 2 < 0, 

2. two active gas with linear isotherms: q*(c\,C2) = KiCi, i= 1,2, 

3. two active gas with binary Langmuir isotherms: ql{c\, C2) = — — — , i — 1,2, where 

1 + K\a + K 2 c 2 

positive constants Q\,Q2, K\ > K 2 satisfy: Q\K\ < Q2^2- 

Furthermore, for two active gas with binary Langmuir isotherms, A is genuinely nonlinear, i. e. h2b)) 
is satisfied, if Q\K\ ^ Q2^2- 

The first case is the most classical case when only one gas is active and his isotherm has no inflexion 
point, for instance the Langmuir isotherm. 

The second case is less interesting in chemistry and only valid when concentrations are near con- 
stant states. 

For the third case, notice that K\ > K 2 is not really an assumption (exchange the labels if neces- 
sary). 

Proof of Proposition I5.lt we use some technical Lemmas postponed to Subsection 15.11 The 
point is to satisfy (f2"9"]l . 

1. Case with an inert gas: we have h = q 2 , f(c) = —ch(c), f — —h ~ ch' , f" = —2h' — ch", 
which implies h' = q 2 < 0, h" = (72" < and then /" > 0. We conclude thanks to Lemmas 15.31 and 

2. Case with linear isotherms: linear isotherms are qi(c) = K\c, 92(c) = i^l — c) with K\ > 
0, K 2 > then q[(c) = Kx > 0, q' 2 {c) = -K 2 < 0, h'(c) = q[(c) + q' 2 {c) = K x - K 2 , f"(c) = 
2(K 2 — Ki). We assume K\ < K 2 , then we have hi < < /" . Since = 0, i = 1, 2, we conclude 
thanks to Lemmas 15.31 and 15.51 



9 n ,, , . , . . , / \ QiKic f . Q 2 K 2 (l-c) 

d. Case with a binary Langmuir isotherm: we have q\[c) — — — — , 92(c) - 



D D 

where D = 1 + K x c + K 2 {1 - c). Then q^c) = > 0, q 2 {c) = — 2 < 0, 
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h'{c) = q[{c) + q' 2 (c) < if and only if Q X K X {1 + K 2 ) < Q 2 K 2 (1 + K x ), 

q{(c) = —5 < if and only if K\> K 2 , 

D 6 

"f ^ 2Q 2 K 2 (1 + K^jK, -K 2 ) 
q 2 (c) = j) 3 > if and only if K\ > K 2 , 

t>«\ 2(Q 2 X 2 - gigi)(l + gi)(l + gg) ^ n , , -fn^^n^ 
/ (c) = > if and only if Q 2 K 2 > Q\K\. 

Since QiK x < Q 2 K 2 , we get /" > and — i < f < \ — ^ because > if 2 , so we have 

Q 2 Ki f + A 2 

< 77 1 , i-e. ft' < 0. Now we conclude with Lemmas 15.31 and 15.51 □ 

Q 2 ~ K x I + K 2 



5.1 Technical lemmas about shock curves 

We express the shock curves as follows. 

Lemma 5.1 We have exp{5(c + , c_)} = — = — — - — , where h± — h(c±) and a = + f . 

u+ a + h+ [c\ 

Proof: first, from the Rankine Hugoniot conditions: - — — \ — j— L i.e. [ft] — L^l£jl ( ? 1 ( c )] we 

[c + qi (c)\ [ft] ' [uc] 

obtain 



«± = [c + gl (c)]- C .[h] (30) 

U- [c+gi( C )]-C+[ft] 

where [c] = c+ — c_ and [ft] = h(c + ) — h(c-) = h + — ft_, and we get (|5D|) thanks to the following 
computations: 

c+g lC -c_ft = c + gl c - c- 1 JL \ WJ = 1 ^ n (H - e-N 

[mcJ [uc] 

= TT7r[ c + 9i(c)], 



r , / m r»ii r , < w M [c + gi(c)] [c + gi(c)] 

i c + <?i(c)J -c+[ft] = [c + gi(c)J-c+ — = — ([uc\~c + [u\) 

"LLC "LLC 



[ C ] U - r / M 



Rewriting ([3H)) we get 



u_ _ [c + ffi(c)] - c + [h] _ [gx] + [c] - c + [ft] _ [gi] + [c] + c + (ft„ - ft + ) 
u + [c + gi(c)] -c_[ft] [gi] + [c] — c_ [ft] [gi] + [c] + c_(ft_ - ft+) 

_ ill] ~ c + /i + + [c] + ft_c + _ [/] + [c] + ft_[c] _ a + ft_ 
[gi] + c_ft_ + [c] — ft + c_ [/] + [c] + ft+ [c] a + h + ' 

which concludes the proof. □ 

dS 

We need to know the sign of a + h± before studying the sign of — — . 

dc± 

Lemma 5.2 If ft' < and c+ < c < c_ fften a + ft(c+) > a + ft(c) > a + ft(c_) > 0. 
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[/] 

Proof: since h < and c+ < c_ we have ft(c+) > ft(c_) and it is enough to show that — + 

1 + ft(c_) > 0. This inequality is equivalent to [/] + [c] + [c]ft(c_) < because [c] = c + — c_ < 0. 
Since /(c) = 91(c) — cft(c) the inequality is equivalent to [qi] + [c] < c + [h}. We know that q[ > 0, 
c + < c_, ft' < then [q x ] < 0, [c] < 0, [ft] > and then [qx] + [c] < < c + [ft]. □ 

Lemma 5.3 if ft' < 0, i// is convex and if c + < c_ then we have — — (c+,c_) < 0. 

oc + 

Proof: we have 5(c+,c_) = [L] = ln(u+) — ln(w_) = ln(-^t) and — — = — — — — - thanks 

u- dc+ u- <9c+ a + ft- 

d cv ~\~ ft 1 ( doc \ 

to Lemma 15.11 A calculus gives — — = — — I — — — [ft] + h'(c + )(a + ft_) I . Now 

oc + a + ft_ {a + h-) z \ oc + J 

dcv 

— — > because / is convex, next [ft] > since h' < and c + < c_ . Lastly a + ft_ > from 
oc+ 

OS 

Lemma [O] and we get — — (c + ,c_) < 0. □ 
oc+ 

The following result concerns the case with an inert gas: 

dS 

Lemma 5.4 If q% =0 and ail < then — — > for c_ > c + . 

OC- 

Proof: if q\ = then /(c) = — cft(c), ft(c) = 52(c) then ft'(c) = q' 2 (c) < 0. By a direct computation 
and thanks to Lemma 1 5. 1[ we have 

u+ _ [c]-c_[ft] _ [c]-c+[h] + [c][h] =1 , 1 



it- [c]-c+[ft] [c]-c+[ft] J__£±' 

[ft] [c] 

9 1 c+ tt+ 

But — — — < 0, —tt decreases, then — increases with respect to c_. □ 
oc_ [ft] [cj w_ 

In the case of two active components we have the following result: 

dS 

Lemma 5.5 If q 1 < < q 2 and if f is convex then — — (c+,c_) > 0. 



Proof: let be c between c+ and c_ . From Lemma I5.2I we have: 

= /(c + )-/(c) + x + > Qj = /(c+)-/(c) + x + > 

c+ — c c+ — c 



We rewrite S using the functions u, v. With Lemma |5. II we get immediately: 

- »m 

The function / is convex, so u is increasing. ^From equality /(c) = qi(c) — cft(c) we have 
(v(c) - l)(c+ - c) = gi(c+) - c + ft(c+) - gi (c) + cft(c) + ft(c)(c+ - c) 
= 9i(c+) -?i(c) -c+(ft(c+) -h(c)). 
Recall that h(c) = qi(c) + 92(c), so we have: 

(v(c) - l)(c+ - c) = 9i(c+) -31(c) -c+(t?i(c+) +9-2(0+) - qi(c) - 92(c)) 
= (! - c +)(<7i( c +) - 9i(c)) - c + (<7 2 (c+) - 92(c)). 

Finally, w(c) - 1 = (1 - c+) gl(c+) ~ gl(c) - C+ ^£±}_^M with < c+ < 1. Now, 9l is concave 

c+ — c c+ — c 

and 92 is convex, so v is decreasing. Finally, — is increasing and — — > 0. □ 

v oc- 
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6 Interactions estimates 

In this section we study the evolution of the total variation of L = ln(u), denoted TVL, through 
waves interactions. It is a key point to obtain some BV bounds and a special structure for velocity. 
Let us denote (cq,L ), {ci,L{), (c2,L 2 ), three constant states such that: 

• the Riemann problem with (cq,Lq) for the left state and (c\,Li) for the right state is solved 
by a simple wave Wi , 

• the Riemann problem with (ci,Li) for the left state and (c 2 ,L 2 ) for the right state is solved 
by a simple wave YV 2 , 

• Wi and W2 interact. 

Just after the interaction we have two outgoing waves W^, VUJ, and the intermediary constant 
state (c*, L\). We denote by TVL the total variation of lnu just before interaction: 

TVL = \Lq — L\\ + \L\ — L 2 \. 
We denote by TVL* the total variation of In it just after the interaction: 

TVL* = \L -Lt\ + \L*-L 2 \. 
We use similar notation for the concentration. 

Denote by a_ the negative part of a: ct_ = max(0, —a) = — min(0, a) > 0. 
We have the following key estimates: 

Theorem 6.1 (Variation on TVlnu and TVc through two waves interaction) 

Assume 126\) . Then there exists T > 0, a true constant such that: 

TVL* < TVL + r|co-ci||ci-c 2 |, (31) 
TVc* < TVc. (32) 

Furthermore, if {27}) is also satisfied then: 

TVL* < TVL + T{cx-cq)- x (ca-ci)-, (33) 

in addition, if S , from if&S)) , satisfies the following triangular inequality: 

S{c 2 , c ) < S(c 2 , ci) + 5(ci, c ) 

when Co > c\ > c 2 , then 

TVL* < TVL. (34) 

Inequality (|32p means that the total variation of c does not increase and Inequality ()33|1 means 
that the total variation of In u does not increase after a wave interaction except when two shocks 
interact. In this last case the increase of TVhiu is quadratic with respect to the concentration 
variation. 

Such estimates are only valid when / has no inflexion point. Else, A-wave curves are only 
Lipschitz and we loose the quadratic control for the total variation of L. 

Proof of Inequality (1 3 2 1) : the decay of the total variation of the concentration is straightfor- 
ward since c is constant through a contact discontinuity i.e. c* = Co : 

TVc* = \c 2 - c\\ + \c\ - c | = |c 2 - c | < \c - ci| + |ci - c 2 | = TVc. □ 

Proof of Inequality ()3ip : this proof is much more complicated. We only assume (f26|) . The 
proof is a consequence of the following lemmas. 

Lemma 6.1 If a X-wave interacts with a contact discontinuity then we have TVL* = TVL. 
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Proof: it is the simplest case. We have c\ = c 2 from the contact discontinuity, so, with T defined 
in (28]), Lx—Lq = T{c\,cq) = T(c 2 , c ) and, since c* = c , we have L 2 - L* = T(c 2 , c*) = T(c 2 , c ). 
Then 

L 2 — L x = L\ — Lo, 

which implies L 2 -L 1 =L\- L and 7V£* = TV"i. □ 

Lemma 6.2 There exists a constant T > suc/i i/iai, /or aZZ Co, ci, c 2 € [0, 1]: 

|T(c 2 ,c ) - T(c 2 ,c x ) - T(d,c )\ |< T | ca -c x || ex - co|. 

Proof: we define i? by J?(a, /3) = T(c 2 , c ) — T(c 2 , Ci) — T(ci, Co). We have to prove that i?(a, /3) = 
0(a(3), where a = c\ — cq, = c\ — c 2 . We denote c = C2, b = ci, a = c . We have T G C 3 ([0, 1],K) 
since A is genuinely nonlinear and T(b, b) = 0. We apply the Taylor's formula: 

T{c,a) = T(b-/3,b + a)=T{b,b)-l3d 1 T(b,b) + ad 2 T{b,b) 

+ f {l-t){/3 2 dfS + a 2 d^T-2af3df 2 T){b~t(3,b + ta)dt, 
Jo 

T(b,a) = T(b,b + a) =T(b,b) + ad 2 T{b 7 b)+ I (1 - t)c?d 2 2 T{b,b + ta)dt, 

Jo 

T(c 7 b) = T(b-0,b)=T(b,b)-pBiT(b,b)+ [ (1 - t)[3 2 d 2 T(b - tp,b)dt, 

Jo 

R(a,f3) = T(c,a)-T(c,b)-T(b,a) 

= -T(b,b)+ f (l-t){p 2 (d 2 T{b~t/3,b + ta)~d 2 T{b~t/3,b)) + 
Jo 

a 2 {dlT{b -t/3,b + to) - b%T(b, b + to)) - 2a/3d 1 d 2 T(b -t/3,b + ta))dt. 

Since 

dlT(b-t/3,b + ta)-dfT(b-t0,b) = 0(ta) = O(o), 
dlT(b-t/3,b + ta)-d%T(b,b + ta) = 0(t(3) = 0(j3), 
d 1 d 2 T(b-tf3,b + ta) = 0(1), 

we conclude that R(a, /3) = 0((3 2 a + a 2 + a/3) = 0(aj3). □ 
To conclude the proof of Inequality (f3Tj) it suffices to use the next lemma. 

Lemma 6.3 If two X-waves interact then we have TVL* < TVL + T\c 2 — ci| \c\ — cq\. 

Proof: by definition of TVL and TVL* it suffices to prove that 

LI = L + 0{\c 2 - Cl \\ Cl - c \), 

since TVL* = \L 2 - L\\ + \L\-L Q \ < \L 2 - L \ + 2\L{ - L \ < TVL + 2\L\ - L \. 

Indeed, we have: L\ — Lq = T(ci, cq), L 2 — L\ = T(c 2 , ci), L 2 — L\ = T{c 2l c*) = T(c 2 , cq). Next: 

L 2 —Lq = T(c 2 , ci) + T(ci, Co) and then 

L\- L = T(c 2 , ci) + T(ci, c ) - T(c 2 , c ), 

which allows us to conclude the proof of Lemma 16.31 with Lemma 16.21 □ 
The proof of Inequality ([3Tj) is now complete. □ 
Proof of Inequalities (|33D, (j34ft : 

we assume again (I27|) and also (1^1) to fix the signs. There are more cases to study: 
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• first, we have yet studied in Lemma 16. II the interaction of a shock wave or a rarefaction wave 
(A- wave) with a contact discontinuity (1-wave): the contact discontinuity is "transparent" 
since TVL* — TVL and the concentration variation is also invariant. 

• second, we study the interaction of a shock wave with a rarefaction wave (A-waves with 
different types): see Lemmas l6.4l 153116. 6l and l6. 71 We get TVL* < TVL and the concentration 
variation decreases. It is the only case where TVL and TVc decrease. 

• finally, we study the interaction of two shock waves. In this situation TVL* > TVL and TVc 
is invariant. 

Furthermore, if S satisfies some "triangular inequality", we get TVL* — TVL. 

In order to simplify the notations we denote by D a contact discontinuity, R a rarefaction wave and 
S a shock wave. " RD — * DR " means that a rarefaction wave coming from the left interacts with 
a contact discontinuity and produces a new left wave, namely a contact discontinuity, and a new 
right wave, namely a rarefaction. 

Since a contact discontinuity has a null speed and a A-wave has a positive speed, the only cases for 

Wi, W 2 are: RD, SD, RS, SR and SS. 

For the resulting waves W*, Wf , there are 7 cases. 

The first two cases RD — > DR and SD — > DS have yet been studied in Lemma [BTTl 
Lemma 6.4 In the case RS -> DR, TVL decreases i.e. TVL* < TVL. 

Proof: at the beginning, we have a rarefaction, then cq < c±, Lq > L%, and a shock, then c 2 < c\, 
L 2 > L\. After the interaction, we have a contact discontinuity, then cq = c*, and a rarefaction, 
then c* < C2, T* > L 2 . Finally, we have cq = c\ < c 2 < Ci then g(co) = g{c*) < g{c 2 ) < g{ci). We 
can write 

TVL = \L -L 1 \ + \L 1 -L 2 \=L -L 1 +L 2 -L U 
TVL* = \L -L* 1 \ + \L 2 -L* l \=\L -L* 1 \+L* 1 -L 2 . 

There are two cases: 

• the simplest is L > T*, then TVL* = L - L\+ L\- L 2 = L - L 2 < L - Ti < TVL, 

• the second case is L < L*. Let us define L 2 by 

Lq — L 2 = L*^ — L 2 , 

then Lq - L 2 = L\ - L 2 = g(c 2 ) - g(c\) = g(c 2 ) - g(c ) < g(c x ) - g(c ) = L - L x because 
[L] = —[g] for a rarefaction and c\ = cq. Since shock curves are decreasing, we know that 
L 2 > L x , so TVL* = L* x -L +L\-L 2 = L 2 -L 2 + L -L 2 < L 2 -L 1 + L -L 1 = TVL. □ 

Lemma 6.5 In the case RS -> DS we get TVL* < TVL. 

dS 

Proof: this case needs the assumption — — > 0. At the beginning, we have a rarefaction: c\ > cq 

oc^ 

and L\ < Lq with a shock: c 2 < c\ and L 2 > L\. The state (c 2 ,L 2 ) is connected with a shock 
(c*, L* ): c 2 < c* and L\ < L 2 . The state (co, To) is connected with a contact discontinuity (c*, T^): 
co = c*. Finally, we have c 2 < cq = c* < c\. Then TVL =| To — Ti | + | L± — L 2 \= Lq — Li+L 2 — Li 

dS 

and TVL* =1 T - L\ I + I T 2 - L\ 1= T, - L*+ \L1-L \. But, with the assumption, - — > 0, 

oc- 

S(c 2 ,cq) = S(c 2 ,cl) = L 2 — T* < S(c 2 , ci) = T 2 — Ti then T* > L\. 
There are two cases: 

• if T > L\ then TVL* = T - L\ + L 2 - L\ < T - Ti + T 2 - Ti = T^T, 

• else T < L\ then TV^T* = -T + L\ + L 2 - L\ = L 2 - L < L 2 - L x < TVL. □ 

Lemma 6.6 In the case SR -> DR we have TVL* < TVL. 



16 



Proof: in the beginning, we have a shock who interacts with a rarefaction then c\ < Co, L\ > Lq 
and c 2 > Cx, Lx > L 2 . 

After the interaction, we have a contact discontinuity then cq = c{ and a rarefaction then C2 > c* 

and > ^2- Finally, we have ci < Co = c* < C2. Since g' > 0, we have g{c\) < g{co) < 5(02). 

For a rarefaction [£] = —[g] then L2 — L\ = g{c\) — g(c 2 ) — g(co) — g{c2) because c\ = cq, 

L 2 — Lx = g(ci) — g(c 2 ) < g{ c o) — g( c 2) because c\ < cq and g' > 0. 

So we have: Li — L\ < g{c\) — gics) = L2 — L\ and 

TVL =| L x - La \ + \ L 2 - L 1 |= Lx - L a + L 1 - L 2 > Lx - L 2 , 

TVL* =| L\ - L I + I L 2 - L\ | = | L\-L Q \ +L\ - L 2 . 

There are two cases: 

• the first is L\ > L then TVL* = L\- L - L 2 + L\ = 2L\-L Q -L 2 = -{L 2 - L\)+L* x - L 
<-(L 2 - Lx) +Ll-L 2 + L 2 -L < -L 2 + Lx - L 2 + L x + L 2 - L = 2L X -L 2 -L a = TVL, 

• the second case is L\ < L then TVL* — —L\ + L — L 2 + LI = L - L 2 < L\ - L 2 < TVL. 

□ 



Lemma 6.7 In the case SR -> DS, TVL decreases i.e. TVL* < TVL. 

This situation is illustrated in Fig. [T] 

Proof: it is the most difficult case. At the beginning, we have a shock then ci < Co and 
Lx > Lq. The shock interacts with a rarefaction then c 2 > cx and L 2 < L\. 
We then have TVL =\ Lx - L | + | L 2 - Lx |= Lx - L + Lx - L 2 . 
The state (C2, L 2 ) is connected to (c*, L\) by a shock then C2 < c\ and L\ < L 2 . 
The state (co,io) is connected to (c\,L\) by a contact discontinuity then cq = c*. 
Finally, we have ci < c 2 < c\ = c , S(cx,c ) = S w > S(c 2 , c ) = = S"(c2, c{) = L 2 — L\, 

dS 

Lx — Lq = Sxo > S 2 q = L-) — LX , because — — < 0. 

oc + 

There are two cases: 

• if L < L* (see Fig. [U left) then L 2 < Lx and 

TVL* =| LJ - L I + I L 2 - LI |= L* 1 -L + L 2 -L* 1 =L 2 -L <Lx-L < TVL, 

• if L* < L (see Fig. O right) then we define L 2 by L 2 — L = S 2 q = L 2 — L* < Sxo = Lx — Lq 
zxidTVL* =\ L\-L I + I L 2 -L\ |= L Q -L\ + L 2 -L\ = L 2 -L 2 + S 20 < Lx-L Q + Lx-L = 
TVL. □ 
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L=ln u 




Figure 2: SR — > DS: first case, left, and second case, right. 



The following case is the only one where TVL increases, except if S satisfies a "triangular 
inequality" . 

Lemma 6.8 In the case SS — > DS we have 

TVL* = TVL + 2max(S 20 ~ S21 - S w ,0) = TVL + 2 max(Lo — L*, 0) > TFL. 

Proof: at the beginning, we have a shock: ci < Co and L\ > Lo- It interacts with an another 
shock: c 2 < C\ and L2 > L\. 

The state (02, L2) is connected to (c*, L*) by a shock then C2 < and L^ < L 2 . 

The state (co, Lq) is connected with with (c*, L*) by a contact discontinuity then cq = c*. 

Finally, we have c 2 < c\ < c — c* and L < L 1 < L 2 - 

With L 2 — Li = S 2 i > 0, L\ — L = S w > 0, L 2 — L\ = S 20 > 0, we have: 

TVL =\L 2 -L 1 \ + \L 1 -L |= L 2 - L x + L x - L = S 21 + S w , 

TVL* =| L 2 - L\ I + I L\ - L | = | S 20 I + I L\ - L 2 + L 2 - L \ = S 20 + \ -S 20 + L 2 - L \ . 
There are two cases to study: 

• if -^20 + L 2 - L > i.e. ^20 = L 2 - L\ < S 2 i + S w = L 2 - L i.e. L < L\ 
then TVL* = S 2Q - S 2Q + L 2 - L = L 2 - L = TVL, 

• else LI < Lo and we have 

TVL = S20 + S20 — L 2 + Lq = 2S 2 q — 2L 2 + 2Lq + L 2 — Lq 

= 2(S 2a - (L 2 - L )) + TVL, = 2(S 2a - S 21 - S w ) + TVL 
= 2(L -Ll) + TVL, 

which conclude the proof of Lemma 16.81 □ 
The proof of Theorem 16. II is now complete. □ 

7 BV estimates with respect to time for the velocity 

In System I©-©-©, there is no partial derivative with respect to t for u. Nevertheless, the 
hyperbolicity of this system ( with x as the evolution variable) suggests that a BV regularity of the 
"initial" data Ub for x — is propagated. Furthermore, in the case with smooth concentration, the 
Riemann invariant uG(c) suggests that when hiub is only in L°°(0,T), we can hope u(t,x)/ub{t) 
to be still BV in time for almost all x. We prove that this BV structure of the velocity is still 
valid with some convexity assumptions, using a Front Tracking Algorithm (FTA). We conjecture 
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that this structure is still valid for the general case without convexity assumption or, better, with 
a piecewise genuinely nonlinear eigenvalue A = H(c)/u. But, in this last case, the FTA becomes 
very complicated (see Dafermos' comments in [20]). 

7.1 The case ln^t 6 BV(0,T) 

Wc first precise the notations used in the next theorem. We define the function cj on (0, T) by 



cj(s) 



c (s) if < s < X 
c b (-s) if < -s < T 



and we set TVa = TVci[-T,X]. 

There exists a positive constant 7 such that if (c_, L_) is connected to (c + , L + ) by a A-wave then 
I L+ — £_ I < 7 I c + — c_ I . That is an easy consequence of (|28|) . Indeed, it is yet proven in [8] , Lemma 
3.1, with an inert gas, or in [S], Lemma 4.1, for two active gases. 
The constant T comes from Theorem 16. II 



Theorem 7.1 [Propagation of BV regularity in time for the velocity] 

Assume (261) . If IrtUb £ BV(0,T), if Co,Q, G BV^ and i/ (u,c) is a weak entropy solution of 
System coming from the Front Tracking Algorithm, then c € BV({0,T) x (0,X)) and 

u e L°°((0,T),BV(0,X)) n L°°((0,X), BV(0,T)). More precisely: 

maxf sup TV x c(t,.)[Q,X], sup TV t c(.,x)[0,T]) < TVa, 

\0<t<T 0<x<X J 

sup TV x \nu(t,.)[0,X] < TV\nu b + jTVa, 

0<t<T 



sup TV t lnu(.,x)[0,T] < TVlnui, + 2jTVci + — (TVci) 2 
o<x<x 2 



Compared to 8, 9j, the new result is that u(t, x) is BV with respect to time if u b is in BV(0, T) i.e. 
the last inequality of the Theorem 17. II With the Godunov scheme used in 8, 9] we do not obtain 
such time regularity for the velocity. It is the reason why we use the FTA to get more precise 
estimates. Notice that we consider a local (in time and space) problem for reasons of realism: we 
could consider a global one as well, i.e. for (t,x) 6 (0, +00) 2 . 

Proof: The easiest BV estimate on the concentration c after interaction (estimate (|32[) in 
Theorem 16. 1[) . which is always valid independently of the velocity u, yields to a control of c in 
LfBV x fl L^BVx as in [S], since A waves always have a positive speed. From Lemma 4.8 of [5] 
p. 80 (ore more simply Lemma 3.1 of [5] p. 557) we get Lf° x fl LfBV x bounds for the velocity 
u. It follows, from a natural adaptation of the estimates and compactness argument of the proof 
of Theorem 5.1 p 563. in [8] or Theorem 6.1 p. 83 in [9], that there exists a subsequence which 
converges to a solution of the initial boundary value problem with the prescribed data co,Cb,Ub 
when S goes to zero, thanks to the approximate entropy inequality (|25|) . Furthermore, as in [3 [9], 
we recover strong traces at t = and x = 0. 

Notice that this existence proof is also valid without any BV assumption on the velocity at the 
boundary: we only need Inub in L°°(0,T). 

The BV estimate with respect to time for In it, i.e. the third estimate in the theorem, is a 
consequence of two following lemmas. □ 

Let (u, c) be an entropy solution coming from FTA. For 8 > 0, representing the distance from 
the boundary x = or t = 0, let us define: 

T , _v _ J \nu(t = \s\,x — 5) if — T < s < 
[Sl > ~ \ \nu{t = 5,x = s) if0<s<V ' 

TVL(0) = limswpTVL(.,S)[-T,X]. 
For piecewise data, TVL(0) is the total variation of lnu just before the first interaction. 
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Lemma 7.1 Before wave-interactions we have TVL(0) < TVlnub + 27 TVc/. 

Proof: it suffices to prove this inequality for a piecewise constant approximate solution issued 
from the FTA. We discretize [0,T] and [Q,X] as follows: 
T = si > s 2 ■ ■ ■ > s m > s m+ i = < s m+2 < ■ ■ ■ < s N — X. 

For i — 1, • • • , m let us define the following piecewise approximations of c and In it: 

1 /" Sl 1 f Sl 

Cb(t)dt, Li — / hx(uf,(t))dt. 



Si Si-|_i J s . + 1 Si Si-|_l 

Since i = is a characteristic boundary we define only Cj for £ = m + 1, • • • ,N— 1 by: 

co(a;)da;. 

Si + l 



For i < m we solve the i" 1 Riemann Problem with left state (q, Li) and right state (cj+i, Lj+i) and 
we denote by c*, L* the intermediary state. Indeed c* = Cj+i since c is constant through a contact 
discontinuity. From Lemma 3.1 p. 557 of [S] (or Lemma 4.1 p. 78-79 of [S] for two active gases) we 
know that: 

\L, -L*\ < 7 | Cl - c*| = 7 | Ci - Cj+i|. 
We now estimate the total variation of lnu for the i th Riemann problem: 

m-Lti + m-L^i < iLi-Lii+im-m+iLi-L^i) 

< 2j\a — Cf+i| + \Li - Li+i\. 

Now, we look at the corner t — 0, x = and i = m. There is only a A-wave since the boundary is 
characteristic. With the left state (c m ,L m ) and only (c m+ i) for the right state, the resolution of 
the Riemann problem gives us a new constant value for In it, namely L m +i = L* n . We have again 
the estimate \L m — L* n \ = \L m — L,„+i| < 7 |c m — c m+ i|. So for i = m + 1, m + 2, • • • ,N — 1 we 
define solving the characteristic Riemann problems with the estimate: 

\Li - L l+1 \ < j\a - Cf+i|. 

Summing up with respect to i, we obtain the total variation on L just before the first wave inter- 
action: 

TVL < 5^(27|c i -c i+ i| + |£ i -Ln. 1 |) + 5^7|ci-c i+ i| 
< TV\nu b + 2-fTVa. 

□ 



Lemma 7.2 We have the following estimate: TVL < TVL(0) + - {TV a) 2 . 

Proof: we prove this estimate for any constant piecewise approximation built from the FTA. The 
same estimate is still true passing to the limit. 

First, we enumerate the absolute value of the jump concentration initial-boundary value from 
the left to the right: 

ai = ^ - Ci_i i = 1, • • • , N. 

Notice that we have N + 1 constant states for the initial-boundary data: (cq, Lq), ■ ■ • , (cjv, £jv)- 

From Theorem 16. 1[ the increase of the total variation of lnu is governed by following inequality 
TVL* < TVL + r|o!j_i I |q!j I if the wave number i — 1 interacts with the wave number i. Since c 
is constant through a contact discontinuity (c is a 2-Riemann invariant) and the jump of c adds 
up if two A-waves interact, we consider only interaction between A-waves. Indeed we neglect that 



20 



interaction with rarefaction has the tendency to reduce TVL. 

We measure the strength of A-wave with the jump of c through the wave. We have positive or 
negative sign whether we have a rarefaction wave or a shock wave. 

Let (pqj 1<i<N _ k be the strength of the A-wave number i (labeled from the left to the right) after 
the interaction number k. We have a? = on and denote by j k the index such that the interaction 
number k occurs with the A-wave number j k and j k + 1 where 1 < jk < N — k. For 1 < i < N — k, 
the strengths of A-waves after the interaction number k > are given by: 

if i < j k 
if i = j k , 
if i > j k 

and the increasing of TVL is less or equal than rS' fc where, from Theorem 16.11 

S° = 0, S k =S k - 1 + \a k - 1 \\a k ^\. 

Let us define the integers l k as follows: 
l\ = i and at each interaction 

, fc f It 1 if i<3 k , 
1 I if • / V /.• • I. 

Notice that after each interactions with two A-waves, there is only one outgoing A-wave. Thus, the 
number of A-waves decreases at each interactions, which proves again (see |21| ) that the number of 
interactions is finite and the FTA is well posed. 

By induction, we see that: a k = a\ where l\ = 1 < 1% < ■ ■ ■ < l%_ k+ i = N — k + 1, 

l ( - = i and l k is non decreasing with respect to k. Now, from the definition of S k , we can deduce 
that: 

S k = S k-i + \ai\\aj\, (35) 

(jj)ej fc 

where J k = {(ij); ifr 1 < i < l k ^ < j < l^}. 
Let us check that: 

s k = Yl HN< ( 36 ) 

(i,j)ei* 

where = 1° C I 1 C • • • C I^ 1 C I k C • • • C / = {(i,j); 1 <i < j < N}. 

It is true for k = 0. It is true for all k if I k ~ l n J k = and then I k = I^ 1 U J k . The point 
is only to prove that I k ~ 1 n J k = 0. Terms la.iHajl in the last sum of ([33]) have indexes i and j 
which appear in two consecutive intervals, i.e. l^ 1 < i < < J < ^fc+ 2 anc ^ ai t er i f° r i = j fc ; 

l k = l k ~ x and l k +1 — l k +2 ■ ^° * an( i 3 ^Ve in the same interval and then terms laiHajl cannot 
appear again in S k+1 , S k+2 , . . . , since such intervals are not decreasing. 

The same is true for all indexes in I k . They can appear at most one time in S k . We then have 
jk-i n Jk = g and q ig true ^ 

We easily estimate S k , which concludes the proof: 

N N 

s k < Yl klK-l<2EEl^lK-l 




<^(T^ C/ ) 2 . 



□ 
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7.2 The case \nu b e L°°(0,T) 

For In Mb G L°° and cq, Cb S we get a £?V structure for the velocity. 
Theorem 7.2 [ BV structure for the velocity] We assume V26}) . 

7/lnuft G L°°(0,T), if Cq, G -BV anci if (c,u) is a weak entropy solution issued from the FTA, 
then 

max) sup TV x c(t, .)[0,X], sup TV t c(., x)[0, T] ] < TVci 

\0<t<T 0<x<X J 

and there exists a function v and constants 7, T > such that u(t,x) = Ub(t) x v(t, x) with 
Inv G {L ao ((0,X),BV(0,T))nL oo ((0,T),BV(0,X))} C BV((0,T) x (0,X)), 

sup TF x lnu(t,.)[0,Jf] < jTV Cl , 

0<t<T 

sup TF t ln«(.,a;)[0,r] < 2 7 TVcj + — (TVci) 2 . 
o<x<x 2 

nit X) 

The new result in this theorem is that }—— is BV with respect to time, although Ub is not 

u b (t) 

assumed to be BV, but just in The other regularity properties have yet been proved in [SUH] • 

Proof: the first estimates for c are easily obtained as in Theorem 17.11 since the total variation 
of the concentration does not increase after an interaction. The existence proof of such entropy 
solution follows the beginning of the proof of Theorem 17.11 which is a natural adaptation of the 
existence proof from [51 [5] with only L°° velocity. 

We now study the new BV estimates for v. We can define v by the relation u(t,x) = Ub(t)v(t, x) 
because Ub > 0. Let be M — Inv and Mb — lnt>(-, x — 0). The initial total variation of M on x = 
is TV M b = since v(t, x = 0) = l. 

We approach Ub with a piecewise constant data (thus in BV) and we show that the BV estimate 
for M is independent of uj,. Notice the fundamental relation: 

[L] — lnu + — lnu_ = ln(iif,(£) v + ) — ln(u(,(t) u_) = lnv + — lnu_ = [Af]. 

The equality [L] = [M] implies that the A- waves (f2"5f are the same in coordinates (c, L) and (c, M). 
Then, Theorem 16.11 is still valid replacing L by M. We then can repeat the proof of Theorem 17.11 
to get BV estimates for v. □ 

8 Weak limit for velocity with BV concentration 

When c is only in BV, we cannot reduce System ((4]) to a scalar conservation law for c as in section 
[3] Indeed, since the shock speeds depend on the velocity, we have a true 2x2 hyperbolic system. 
Nevertheless we can state following stability result. 

Theorem 8.1 (Stability with respect to weak limit for the velocity ) 

Let (hr(it^)) 0<£<1 be a bounded sequence in L°°(0,T), such that 

u% u b in L°°(0,T) weak *. 

Let be cq G SV((0, X), [0, 1]) and Cb G BV((0, T), [0, 1]). Let (c e , u £ ) be a weak entropy solution of 
System Op on (0, T) x (0, X) issuing from the FTA with initial and boundary values: 

c £ (0,x) = c (x), X > x > 0, 
< c e (i,0) = Cb(i), T>t>0, 
u £ (t,0) = u £ b {t), T>t>0. 
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Then, there exists (u(t,x),c{t,x)), weak entropy solution of System ^ supplemented by initial and 
boundary values: 

c(0, x) = cq(x), x > 0, 
c(t,0) = c b (t), t>0, 
u(t,0) = u b {t), t>0, 
such that, when e goes to and up to a subsequence: 

c e {t, x) -> c(t, x) strongly in L x ([0, T] x [0, X]), 
u e (i,x) u{t,x) weakly in L°°([0,T] x [0,X]) weak *, 

u £ (t,x) = x u(t,a;) +o(l) strongly in L 1 ([0,T] x [0,X]), w/iere u(i,a:) = ^ 



u b (t) 



For the convergence of the whole sequence we need the uniqueness of the entropy solution for 
initial-boundary value problem: (01, ((3J). 

Proof: from Theorem 1 7. 2 1 we know that u £ (t, x) = u\(t)v £ (t, x) where the sequences (lni> e )o< e and 
(c £ )o< e are uniformly bounded in BV((0,T) x (0, X)). Then, up to a subsequence, we have the 
following strong convergence in L 1 ((0, T) x (0, X)): v £ — > v, c £ — > c. 

(c £ ,u £ ) is a weak entropy solution for (j4|) means for all rp such that '0" > and Q such that 
Q' = h't/j + Hijj' we have in distribution sense: d x (u £ (t,x)ip(c £ )) + dtQ(c £ ) < 0, which is rewritten 
as follows: d x (ul(t)v £ (t, x)i/j(c £ )) + dtQ(c £ ) < 0. Passing again to the weak-limit against a strong 
limit we get: d x (u b (t)v(t, x)ip(c)) + dtQ(c) < 0. i.e. (c, u = u b x v) is a weak entropy solution for 
System (j4]). We also can pass to the limit on initial-boundary data. 

Since there exists 8 such that < 5 < u £ b < S^ 1 , v £ (t,x) — > v(t,x) means u £ (t,x)/ul(t) — 
v(t, x) — + and also means u £ (t, x) — u\{t) x v(t, x) — » 0, which concludes the proof. □ 

An example of high oscillations for velocity: as an example of weak limit we consider the case 



of high oscillations for velocity on the boundary. 
Let be u b {t,9) G L°°((0, T), C°(R/Z, R)), u 6 (t) = 



/ u b (t,8)d0 and assume infuj > 0. With 



ul(t) = u b ft, - J and the same notations as in Theorem 18. II we have: 

• first, oscillations do not affect the behavior of the concentration since (c e ) converges strongly 
in L 1 towards c and the limiting system depends only on the average u b and not on oscillations; 

• second, (u £ ) converges weakly towards u b {t) x v(t, x) and we have a strong profile for u £ : 

= 0, where U(t, x, 6) = u b (t, 6) x v(t, x). 



lim 



u £ (t, x) — U I t, x, - 



s 



i 1 ((0,T)x(0,X)) 
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